Search at the CERN LHC for a light neutralino of cosmological interest Q 
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We address the problem of a search at the LHC for a neutralino whose mass is around 10 GeV, 
i.e. in the range of interest for present data of direct search for dark matter particles in the galactic 
halo. This light neutralino is here implemented in an effective Minimal Supersymmetric extension 
of the Standard Model at the electroweak scale without requirement of a gaugino-mass unification 
at a grand unification scale. Within this model we select a representative benchmark and determine 
its prospects of reconstructing the main features of the model at different stages of the LHC runs. 

PACS numbers: 95.35.+d,11.30.Pb,12.60. Jv,95.30.Cq 



^ ■ I. INTRODUCTION 

Most investigations on the search for neutralinos at the CERN Large Hadron Collider (LHC) concern neutralinos 
within Supcrgravity-inspired (SUGRA) models. This implies, in particular, neutralinos of a mass larger than about 
50 GeV, since this is the lower bound on the neutralino mass which directly follows from the LEP lower bound on 
the chargino mass combined with the assumption, inherent in SUGRA models, that the gaugino masses are unified 
at a grand unification (GUT) scale. 

Relaxation of this hypothesis of gaugino-mass unification allows the neutralino mass m x to be smaller than 50 GeV. 
A supersymmetric model which incorporates this possibility and is still very manageable, since expressible in terms 
of a limited number of independent parameter, is the one considered in Ref. It consists of an effective Minimal 
Supersymmetric extension of the Standard Model (MSSM) at the electroweak (EW) scale, whose main properties are 



0^ \ summarized in Sect. [II] 

In Ref. [lj it was stressed that, in case of R-parity conservation, a light neutralino (i.e. a neutralino with m x < 50 
GeV), when happens to be the Lightest Supersymmetric Particle (LSP), constitutes an extremely interesting candidate 
for the dark matter (DM) in the Universe, with direct detection rates accessible to experiments of present generation. 
In Ref. [l| it was also derived a lower bound for m x from the cosmological upper limit on the cold dark matter 
?\ ' density; the value of this lower bound, updated on the basis of the experimental data available in Autumn 2010, was 
established in Ref. [3| to be m x ~ 7.5 GeV. The implementation of the very recent upper bound on the branching 
ratio for the process BR(B® — > (i~ + /i + ) 01 moves now the lower bound on the neutralino mass to the value m x ~ 9 
GeV. The theoretical framework which allows neutralinos with a mass in the range 9 GeV < m x < 50 GeV is briefly 
summarized in Sect. II and will be simply denominated as Light Neutralino Model (LNM) in the present paper. 

In Ref. 0| the LNM was proved to fit the annual modulation effect measured by the DAMA collaboration |5| . Our 
model is also compatible with all experimental searches for indirect evidence of SUSY and with precision data that 
set constraints on possible effects due to supersymmetry, as discussed in detail in Ref. Q (for the compatibility of very 
light neutralino masses with various laboratory bounds see also Ref. Q). Moreover, the possible impact of some early 
analyses by the CMS and ATLAS Collaborations at the LHC on the LNM was investigated in Ref. ■ 

At the same time much interest has recently been raised by a new measurement of an annual modulation effect 
by the CoGcNT Collaboration Q and by some hints of possible signals of dark matter (DM) particles in other 
experiments of direct detection (CDMS 0, CRESST [loj]). 
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What is intriguing in all the experimental results listed above is that, if actually due to a Weakly Interacting 
Massive Particle (WIMP) with a coherent interaction with the atomic nuclei of the detector material, they would all 
be explained by a WIMP physical region with a light mass (around 10 GeV) and a nucleon elastic cross-section in 
agreement with the intervals for the same parameters established by the DAMA Collaboration from a measurement of 
the annual modulation that has now reached a high statistical significance by a running over 13 yearly cycles with the 



DAMA/Nal and the DAMA/LIBRA experiments [ll|. Specifically, compatibility of the DAMA result with CDMS 
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has been discussed in Ref. [12j . and between CoGeNT and DAMA in Ref. [13l |. 

These results have prompted a large number of phenomenological papers focused on WIMPs with a light mass 
Turning to a specific candidate, it has now become common to consider ncutralinos with a mass of order 10 GeV. 

Due to the relevant role that light ncutralinos can have in cosmology and astrophysics it becomes so of the upmost 
interest to investigate the possibility of searching for these particles at the LHC A preliminary analysis in this 



direction was performed in Ref. 15j : event rates were determined in specific scenarios and benchmarks dictated by 
the relevant cosmological properties of the LNM. There it was shown that the perspectives of a fruitful investigation 
of the supersymmctric parameter space relevant for light ncutralinos at the LHC are potentially good, though no 
specific analyses of the signal/background ratios and of kinematical distributions were performed. 



In the present paper we wish to extend the investigation of Ref. |15| by making use of a numerical simulation to 
estimate in a realistic way the detectability of the LNM at the LHC over the SM background and to show what 
information about the masses of SUSY particles can be extracted from the data. 

The paper is organized as follows. In Section |H] we summarize the main features of the LNM, concentrating on the 
particle mass spectra relevant for our analysis. In Section Mil we discuss the general properties in the LNM of the 
decay chains that are used to reconstruct the SUSY masses. In Section HVl we explore the LNM parameter space to 
select a suitable benchmark for our scenario, that is then used in Section[V]to assess its chances for an early discovery 
within the 7 TeV run at the LHC, and in Section IVII to discuss how the SUSY masses can be reconstructed from 
invariant-mass spectra measured with a larger collected luminosity at 14 TeV. We then give our conclusions in Section 

Em 



II. A MODEL FOR LIGHT NEUTRALINOS (LNM) 



The supersymmetric scheme we employ in the present paper is an effective MSSM scheme (effMSSM) at the 
electroweak scale, with the following independent parameters: Mi, M2, M3, /x, tan/3, vtla, mq, mi, mj and A. Notations 
are as follows: Mi, M2 and M3 are the U(l), SU(2) and SU(3) gaugino masses (these parameters are taken here to be 
positive), fi is the Higgs mixing mass parameter, tan/3 the ratio of the two Higgs v.e.v.'s, ttia the mass of the CP-odd 
neutral Higgs boson, is a squark soft-mass common to the squarks of the first two families, mi is the squark 
soft-mass of the third family, mj is a slepton soft-mass common to all sleptons, and A is a common dimensionless 
trilincar parameter for the third family, = A^ = Ami an d -^f = Ami (^ ne trilinear parameters for the other families 
being set equal to zero) . In our model, no gaugino mass unification at a Grand Unified (GUT) scale is assumed: this 
implies that Mi and M2 are independent parameters at the EW scale. The model introduced here is the one discussed 
in Ref. [l| , with the minimal extension that the degeneracy between the soft squark mass of the first two families and 
that of the third family is removed. In particular, the splitting between mq and reduces some tuning introduced 
in the parameter space by the constraint on the b — > s + 7 when mi=niq > 700 GeV Q . 

The following experimental constraints are imposed: accelerators data on su pers ymmetric and Higgs boson searches 



(CERN e + e collider LEP2 |16(, Collider Detectors DO and CDF at Fermilab (171); early bounds from Higgs searches 



at the LHC Ha, LL9J; measurements of the b — >• s + 7 decay process [20|: 2.89 < B{b — > s + 7) • 10 4 < 4.21 is employed 
here (this interval is larger by 25% with respect to the experimental determination [2p| in order to take into account 
theoretical uncertainties in the supersymmetric (SUSY) contributions [2l[ to the branching ratio of the process (for 



the SM calculation, we employ the recent NNLO results from Ref. |22|)); the measurements of the muon anomalous 

of the experimental world average from 



7 cxp 



,the 



magnetic moment a M = (g^ — 2)/2: for the deviation, Aa^ 
the theoretical evaluation within the SM we use here the (2 a) range 31 < Aa M • 10 11 < 479, derived from the 



latest experimental [231 ] and theoretical 24j data (the supersymmetric contributions to the muon anomalous magnetic 
moment within the MSSM are evaluated here by using the formulae in Ref. |25[); the upper bound on the branching 
ratio BR{B° S ->• uT + /z+) [1 0]: we take BR(B° S -> fi~ + < 1.5 ■ 10~ 8 ; the constraints related to AM B , S = 
M B , - M Ss the measurements of the decays B -> tv Q and R(D) = BR(B -> Dtv)/BR(B -> Dev) (30 1 
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(in particular, the compatibility of very light ncutralino masses with the latter four constraints is discussed in detail 
in Ref. 0). 

The linear superpositions of bino B, wino and of the two Higgsino states ff 2 which define the four 

neutralino states, Xi (i = 1j 2, 3, 4) are written here as: 



V 



= o^B + afw^ + 4 } H? + af # 2 °. (1) 



The properties of these states have been investigated in detail, analytically and numerically, in Ref. |15[ for the case 
when the smallest mass eigenstate xi (or x i n short) is light, i.e. m x = m Xl < 50 GeV. Of that analysis we report 
here only the main points which are relevant for the present paper. 



We first notice that the lowest value for m x occurs when: 



m x ~M X « \lA,M 2 . (2) 

since the LEP lower limit on the chargino mass (m x ± > 100 GeV) sets a lower bound on both and M 2 : |/x|,M 2 > 
100 GeV, whereas M\ is unbound. Thus, x = Xi is rnainly a Bino, whose mixings with the other interaction eigenstatcs 
are given by: 



/ (i> 
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(3) 



MiS/3 + /zcg ' 



where £i = mi — M\. These expressions readily follow from the general analytical formulae given in Ref. [1 51 ] by 
taking tan/? > 10, as consistent with the scenarios discussed below. 

Useful approximate expressions obtain also for the compositions of the eigenstates corresponding to the asymptotic 
mass eigenvalues: m; ~ ±/i and rm ~ M 2 . That is: 

a) for the neutralino states Xi with rm — ±M; 



-cot e , 



M2T/1 

2g 2 s e (± M -M 2 ) 
M zSp {se 2 M 2 T^y 



(4) 



where £2 = ±A* — m i! 

b) for the neutralino state Xi with rrij 



M 2 



— tang, 
1VI2 

&a (M| - M 2 ) 
M z {M 2 cp + [isp)c 2 e M 2 
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(5) 
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where £3 = M 2 — rrii . 

From the above expressions the following relevant properties hold: (i) xi i g mainly a B-ino whose mixing with H° 
is sizable at small /i, (ii) \3 nas a mass |ttt,3 | ~ with a large — iJ| mixing, independently of M 2 , (hi) xi an d 
Xa interchange their main structures depending on the value of the ratio |//|/Af 2 : X2 is dominantly a W-ino (with 
a sizable subdominance of H^) for Mi « and a maximal H\ — if| admixture for M 2 >> whereas X4 is a 
maximal -ff^ — £f 2 admixture for Af 2 << and a very pure W-ino for M 2 >> |/x|. 

For specific spectroscopic schemes, characterized by various internal hierarchies, we will use the denominations 



already introduced in Ref. [151 ] . i.e.: (i) normal hierarchical scheme when M 2 < (ii) degenerate scheme when 



M2 ~ (hi) inverted hierarchical scheme when M 2 > (notice that we always assume (Mi << M 2 , |^| 



A. Cosmologically inspired scenarios 

If light ncutralinos are present in the Universe as relic particles, their abundance Vt x h 2 has to be smaller than 
the observed upper bound for cold dark matter (CDM), i.e. Q x h? < (SlcDMh 2 ) m ax = 0.122 (this numerical value 



represents the 2a upper bound to (^cDAih ) ma x derived from the results of Ref. |31|)- 

This requirement implies a lower limit on the neutralino pair annihilation cross section (T ann through the usual 
expression: 

g*(Xf) 1 < OannV > 

where < a ann v > = Xf(a ann u)i n t, (oWn v)i nt being the integral from the present temperature up to the freeze-out 
temperature Tf of the thermally averaged product of the annihilation cross-section times the relative velocity of 
a pair of neutralinos, Xf is defined as Xf = m x /Tf and g+(xf) denotes the relativistic degrees of freedom of the 
thermodynamic bath at Xf. 

The lower bound on cr ann , implied by the cosmological upper limit on CDM, combined with the constraints due 
to accelerator data and other precision measurements, restricts markedly the overall size of the supcrsymmetric 
parameter space, as depicted in Figs. 1-2 of the second paper of Ref. [l|. In particular, it is instrumental in placing 
the fore mentioned limit m x ~ 9 GeV. 

We thus arrive at the formulation of two specific physical scenarios for the case of light ncutralinos of cosmological 



interest, as delineated in Ref. [15|. These scenarios are determined by the different ranges of the mass vha of the 
CP-odd neutral Higgs boson, and are summarized in Table Q] 

We have a Scenario A when 90 GeV < m A < (200-300) GeV (we recall that m A >90 GeV is the LEP lower 
bound). When ttia is in this range, then the neutralino mass can be as small as ~ 9 GeV, since the cosmological 
upper bound is satisfied due to a sizable contribution to the neutralino pair annihilation cross section by the exchange 
of the A Higgs boson in the s channel. For this to be so, the B-ino component of the Xi configuration must be 
maximally mixed with the H° component (i. e. a^ 1 /a>i — 0.4). From the second expression in Eq. Q one sees 
that this condition is satisfied when /i is small (|/i| ~ 100-200 GeV). Moreover, it turns out that tan/3 must be large 
(tan/3 ~ 30 - 45). The trilincar coupling is only mildly constrained to stay in the interval —1 < A < +1; the slcpton 
soft mass mr and the squark soft mass rriq are unconstrained. In this scenario, the following hierarchy holds for the 
coefficients o^ 1 ' of xi- 



|4 1) i>|a| 1) |»|4 1) i,iai 1) i ) (7) 

as easily derivable from Eqs. (j4j. 

When rriA > (200-300) GeV, the cosmological lower bound on cr ann can be satisfied by a pair annihilation process 
which proceeds through an efficient stau-exchange contribution (in the t, u channels). This requires that: (i) the stau 
mass rrif is sufficiently light, nif ^ 90 GeV (notice that the current experimental limit is rrif ~ 87 GeV) and (ii) xi 
is a very pure B-ino (i.e. (1 — aij 1 ) = O(10~ 3 ). If this is the case, then light neutralinos can exist, but with a mass 
above ~ 15-18 GeV [l|,[32|. As discussed in Ref. [IB], conditions (i) and (ii) require that \fi\ > 500 GeV, tan/3 < 
20; m ; ~ > (100 - 200) GeV; the parameter A is typically in the range —2.5 < A < +2.5, the other supcrsymmetric 
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scenario 


Mi [GeV] 


M [GeV] 


tan /3 


m A [GeV] 


m,- [GeV] 


A 


~ 10-14 


110-140 


30-45 


~ 90-110 




B 


~ 25 


> 500 


< 20 


> 200 


100-200 



TABLE I: Representative features for scenarios A and B described in Section Til Al In scenario A: -1 < A < +1, in scenario B: 
-2 < A < +2. 
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FIG. 1: Topologies of the decay q — > qyll. Left: sequential decay. Right: branched decay. 



parameters are not a priori fixed. The sector of the supersymmetric parameter space characterized by these features 
is denoted as Scenario B. Within this scenario it follows from Eqs. (H|) that the following hierarchy holds for the 
coefficients of xi- 



>> a 



ID 



3 h I "2 M"4 



(8) 



Table I summarizes the representative features of scenarios A and B to be used below for the definition of our 
benchmark. 

In the present paper we focus our investigation of scenario A, postponing the discussion of scenario B for a subsequent 
publication. 

III. SIGNALS AT THE LHC 

Squarks and gluinos are expected to be copiously produced in the pp scattering processes at the LHC: pp — > 
qq, qq* ,gg, qg. In turn, squarks, produced cither directly or through gluinos, can generate the sequential decay chains: 



q^ q\t -> qff -> qffxi, 



(9) 



and the branched ones: 



q -> QXi -> q( z , K H, A)xx -> qffxi, 



(10) 



where / stands for a fermion, / for its supersymmetric partner; from now on the neutralino subscript i can only take 
the values 2, 3 or 4. These two topologies are shown in FigfT] 

These are the key processes to be studied at the LHC to measure the sfcrmions and ncutralinos masses jaUlil. They 
would be characterized by hard jets, specific two-body decays and a transverse missing energy (under the hypothesis 
of R-parity conservation). The determination of the masses cannot proceed through a full reconstruction of the decay 
chains, since the LSP neutralino escapes detection, but rather by measurements of specific features in unidimensional 
and multidimensional distributions in kincmatical variables. 

Typical strategics for determining the sfermion and neutralino masses consist in: a) measurements of endpoints in 
single invariant mass distributions 33-38J, b) correlations among different invariant mass distributions 



A. General properties of the decay chains 

The decay chains (|9"))- (fl"U]) have in common the first step, i.e. the squark decay q — > qxi, which can proceed either 
through gauge coupling (which involve the gaugino components of Xi), or Yukawa coupling (which involve the higgsino 



6 



components of Xi)- In the following we will assume for simplicity a situation where the gluino is decoupled, by taking 
M 3 ^> rriq. In this case squarks can be only produced with the same flavor of the partons inside the protons which 
induce the hadronic processes at the LHC. As a consequence of this Yukawa couplings have a subdominant role as 
compared to the gauge couplings, since the relative importance of the Yukawa couplings to the gauge ones depends on 
the ratio m q /mz (jn q and mz being the quark mass and the Z-boson mass, respectively) and heavy flavors are scarce 
in the proton composition. In particular, in this case in the process q —¥ qxi, the Xi' s having a dominant gaugino 
composition are preferentially produced. 

Sequential chains are differentiated from the branched chains by the features of the decay process undertaken by the 
intermediate neutralino state Xi- I n the sequential chain the decay proceeds through the process: Xi ff f fXi 
with a branching ratio BR(xi —> ff —> ffxi) = BR(xi //) BR(/ — ^ fxi)- I 11 the following we will limit our 
considerations to the most interesting cases, where / is a charged lepton (i.e. f = I = e, fJb,r). The size of BR(xi — > 11) 
depends sensitively on the Xi composition. If Xi is dominantly a gaugino, because of the universality of the gaugino 
couplings, the branching ratios BR(x% — ► 11) for the three lepton flavours are about the same ; if Xi is dominantly a 
Higgsino, Xi decays predominantly into a f r pair. 

In the branched chain, Xi decays either through the Z-boson or through a Higgs boson. The first case, i. e. 
Xi — > Z + xi, involves only the Higgsino components of the two neutralino states; the Z boson subsequently decays 
into all (kinematically possible) // pairs according to the Standard Model branching fractions. The second case, i. 
e. Xi {h, A, H) +Xij i n order to have a sizable BR, requires that one neutralino state is dominantly a gaugino, the 
other dominantly a Higgsino. Since in the scenarios considered in the present paper xi is dominantly a B-ino state, 
Xi — > (h, A, H) + xi is of interest when Xi is dominated by the Higgsino components. Because of the hierarchical 
character of the Yukawa coupling, the subsequent decays of the Higgs bosons are dominated by the production of a & 
- b pair. 

A detailed discussion of the branching ratios for the various processes involved in the decay chains (j§j)- (fTU|) for the 
LNM are given in Rcf. fl5| . 



IV. SELECTING A BENCHMARK IN THE LNM SCENARIO 



In the present paper we wish to address the following two points: i) Is the LNM scenario A (as defined in Table Q} 
detectable at the LHC over the SM background? ii) Is it possible to kinematically reconstruct the neutralino mass at 
the LHC in its low range within the LNM scenario (m x ~ 10 GeV)? 

As discussed in Section III A[ one of the basic features of the LNM scenarios is that some of the SUSY parameters 
are forced into rather strict intervals. As shown in Table HI this is particularly true in Scenario A for the parameters 
Mi, /i, tan j3 and tua, the latter three parameters being just beyond the LEP and Tevatron sensitivities. 

A scatter plot of these 4 parameters where all the experimental constraints listed in Section |H] are implemented 
is given in Fig|2l here the parameters have been varied in the following narrow ranges: 10 GeV < M\ < 14 GeV, 
110 GeV < fi < 140 GeV, 30 < tan/3 < 40, 90 GeV < m A < 105 GeV, in accord with the intervals of scenario A in 
Table HI As far as these parameters are concerned, the choice of a benchmark is quite restricted. 

In Fig. [2] we plot with a triangular symbol the representative point that we adopt in Table |TT] as our benchmark: 
Mi=14 GeV (corresponding to m x ~ 11 GeV), fi— 126 GcV,tan/3=34 and toa=97 GeV. The LNM scenario is basically 
independent on the remaining six parameters of the model (M2, M3, rriq, mj, mj and A), which are only constrained 
by the various experimental limits listed in Section |H] In particular LHC physics is very sensitive through the SUSY 
production cross section to the rriq parameter (which drives the mass of squarks of the first two families corresponding 
to the flavors more abundant in colliding protons) and to the gluino mass M3 . The LHC early runs have already started 



to introduce constraints on these parameters, which, however, strongly depend on the adopted SUSY scenario [40h44 [ . 
Following the approach already adopted in Ref. (l5| , for simplicity in the present paper we will limit our discussion to 
the case in which the gluino is heavier than the squark, and for definiteness we will fix it at the representative value 
M 3 = 2 TeV. 

In Ref. the possible impact of some early analyses by the CMS and ATLAS Collaborations at the LHC 40 - 42 on 
the LNM scenario was investigated. The data considered there consisted in the results of searches for supersymm etry 
in proton-proton collisions at a center-of-mass energy of 7 TeV with an integrated luminosity of 35 pb _1 [40j, 
i.e. the results of the CMS Collaboration for events with jets and missing transverse energy [40J, and those of the 
ATLAS Collaboration by studying final states containing jets, missing transverse energy, either with an isolated 
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Mj (GeV) mi (GeV) 

FIG. 2: Scatter plot of the parameters M\ and /j, (left) and of the parameters rriA and tan (3 (right) in the LNM scenario. The 
dots, in black and in gray, correspond to the range of parameters given in the text in accord with scenario A of Table [I] the 
subset of dots in gray show the configurations which fall into region I of Fig. [4] allowing sequential decays through production 
and decay of a X4 (see text). The (red) triangle shows the LNM-seq benchmark given in Table HT1 that is analyzed in detail in 
the present paper. In the right-hand plot the yellow shaded area shows the region disallowed in the plane (m^-tan/?) from 
the results of Refs. as derived in the analysis of Ref. 
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FIG. 3: Squark production cross section at the LHC as a function of the mass of the squarks of the first two families 
m aqua rk,i2 — Triq for a center-of-mass energy Ecm~7 TeV and for a heavy gluino mass, M$=2 TeV. The solid line shows the 
next-to-leading order value calculated with Prospino [45| while the dashed line shows the same quantity at the leading order. 
The highest horizontal line marked "35 pb _1 " shows the bound on the production cross section from early analyses of CMS 
and ATLAS Collaborations at the LHC as calculated in Ref. [?| for the LNM scenario. The lower horizontal lines show an 
estimation for the sensitivity of the LHC run at Ecu =7 TeV for two different values of the collected luminosity. The vertical 
line at m aquar k=7GQ GeV indicates the value of the parameter rriq ~ m squar k,i2 that is adopted in the LNM-seq benchmark 
introduced in Table ILT1 
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FIG. 4: Scatter plot of the parameters Mi and mf in the LNM scenario with m«j=700 GeV and Ms=2 TeV. All plotted 
configurations are allowed by the observational constraints summarized in Section [TT] In this plot the parameters Mi, /i, tan/3 
and m,A are sampled in the ranges of Scenario A indicated in Table [I] and their numerical scan is shown in Fig. [2] The 
remaining two parameters of the model, A and m t -, are shown in Fig(5] All the configurations are subject to the experimental 
bounds listed in Section (TT] The points are plotted in 5 increasingly dark tones of gray corresponding to the following sub- 
intervals for the allowed range for the muon anomalous mag netic moment: 31< Aa M x 10 11 <130, 130< Aa M x 10 11 <230, 
23CK Aa M x 10 11 <330, 330< Aa M x 10 11 <430 and 430< Aa M x 10 11 <479. The regions indicated by numbers I, II, III, IV, 
V and VI correspond to different kinematic regimes for the sequential decay (see text for details). In particular in regions I, II 
III and IV sequential decays are kinematically accessible through a next-to-lightest neutralino \i with i=4 (region I), i=3,4 
(region II), i=2,3,4 (region III) and i=2,3 (region IV). The triangular symbol indicates the LNM-seq benchmark given in Table 
[IT] that is analyzed in detail in the following Sections. 



lepton (electron or muon) 41 or without final leptons 42 1. As reported in Refs. 40, 41 1 the data appeared to be 
consistent with the expected Standard Model (SM) backgrounds; thus an upper bound on the SUSY production cross 
section er(pp — > squark) at the LHC for a center-of-mass energy Ecm=7 TeV was derived in Ref. [7| for the LNM 
scenario. In Fig. [3] we plot this cross section as a function of the common mass of the squark of the first two families 
fnsquarki2 — mq for M%=2 TeV. Moreover, the highest horizontal line marked "35 pb _1 " shows the upper bound on 
the same quantity as derived in Ref. Q, that implies a lower limit m,q > 450 GeV when the gluino is heavy. In the 
same figure the lower horizontal lines show an estimation for the sensitivity of the LHC run at Ecm=7 TeV for two 
different values of the collected luminosity, C = 1 fb _1 and C = 5 ft) -1 , naively obtained by scaling down the bound 
on the cross section from [t| with the square root of the exposition. Assuming C = 5 fb" 1 as the expected total 
collected exposition in the LHC run at Ecm~7 TeV before the stop scheduled for the end of the year 2011, one can 
see that a value rriq ~700 GeV would be by that time just on the verge of discovery, possibly already providing a 
small excess over the standard model background. For this reason in the following we will adopt mq=700 GeV as our 
benchmark value of the soft squark parameter for the first two families. 

We proceed now to discuss the remaining 4 parameters, M2,nij, m,f and A. As discussed in Section Till Al the 
properties of the sequential decay of Eq. ([9]) that we wish to analyze depend sensitively on the hierarchy among the 
masses of the particles involved, i.e. on the masses of the squarks, of the next-to-lightest neutralino and of the slepton. 
In our scenario these three mass scales are driven by m g -, mi and M% or determining in particular whether the 
spectrum of neutralinos is normal or inverted and if the decay Xi=2,3,4 — ► U is kinematically allowed. Since we have 
fixed m,j=700 GeV and in Scenario A the /i parameter is constrained to the narrow range 110 GeV< [i < 140 GeV, 
in Fig. Uwe discuss this mass hierarchy in the plane of the remaining two parameters, M2 and mj. In this plane we 
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schematically represent with a line at the constant value mr=140 GeV the scale of the /j parameter. Then, with the 
exception of the narrow band where M 2 ~ /i, one has fi ~ rn X2 a , M 2 ~ X4 when M 2 > /x (normal hierarchy) and 
(i ~ m X3 4 , M 2 ~ m X2 when M 2 > mu (inverted hierarchy). Moreover mr fixes the scale of the slepton masses. This 
implies that schematically one can divide the M 2 -mr plane in six regions: 

• Region I (M 2 > mr > fi). Here only the decay xi ~> lS kinematically allowed, with m X4 ~ M 2 . 

• Region II (/i > mr > M 2 ). Here only the decays X3A ~~ ^ U are kinematically allowed, with m xa 4 — /i; 

• Region III (M 2 > fi > mr). All decays X2,3,4 — >■ ^ are kinematically allowed, 

• Region IV (M 2 > m,j=700 GeV, /i > mi). Here only the decays X2,3 —> U arc kinematically allowed, with 

m X2,3 — Mi 

• Regions V and VI. No sequential decays are kinematically allowed. 

In the same figure the scatter plot represents a scan of the LNM parameter space with 110 GeV< M 2 < 1500 GeV, 
110 GeV < m? < 1500 GeV, while the other parameters are in the ranges of Scenario A given in Table [U The points 
are plotted in 5 increasingly dark tones of gray corresponding to the following sub-intervals for the allowed range 
for the muon anomalous magnetic moment: 31< Aa p x 10 11 <130, 130< Aa M x 10 11 <230, 230< Aa p x 10 11 <330, 
330< Aa M x 10 11 <430 and 430< Aa M x 10 11 <480. In this way it is possible to see that the kinematic regions 
II and III (at least for the particular choice of m$ adopted here) are not allowed by the upper bound on Aa^ 1 . 
This restricts the present discussion to the possibility of having a sequential decays only in regions I and IV. In 
region IV sequential decays proceed through production and decays of very light next-to-lightest neutralinos X2,3 with 
mr < m Xl2 ~ [i < 140 GeV that, according to the discussion of Section [Hi arc of higgsino type. These features make 
the detection of sequential decays in this case quite challenging, since the corresponding branching ratio is suppressed 
both by the small available phase space, and by the fact that, due to their higgsino nature, a large fraction of the X2,3 
particles decay through the branched topology to a Z or a Higgs boson (see Fig. [TJ right). Moreover, the higgsino 
nature of the X2.3 particles also implies that when a sequential decay actually takes place it mainly proceeds to tau 
final states that are more difficult to measure compared to muons and electrons. As a consequence of this, for the 
choice of a benchmark for our discussion of sequential decays we decide to focus on the only remaining possibility, i.e. 
region I. In Figs. [2] and [5] the configurations belonging to region 1 are plotted in grey. 

Region I looks more promising than region IV for sequential decays. In fact in this case decays proceed through a Xi 
with m Xi ~ M 2 which is of Wino type implying a smaller coupling to the Z and Higgs bosons which reduces branched 
decays and leading to comparable signals to electrons, muons and taus in sequential ones. Moreover, depending on the 
choice of the M 2 and mj parameters, the phase space available to both decays q — > qxi and Xi ~~ ^ ^1 ca n be sizeable. 
As the scatter plot of Figure [4] shows, this still allows for a wide range of possibilities. With the spirit of choosing 
light values for both the Xi and the sleptons masses and to maximize at the same time the phase space available to 
the decay, in Table [TTI we choose as our benchmark the values M 2 =500 GeV and 771^=300 GeV, corresponding to the 
point shown in Fig. [4] with a triangular symbol lying somewhat in the center of the triangle of region I and close to 
the lower bound on both parameters from the muon anomalous magnetic moment. 

We conclude the present discussion with the last two parameters, A and mz. Apparently, since they do not affect 
directly either the neutralino relic abundance in the LNM scenario or the sequential decays we wish to discuss, pinning 
down a value for these two parameters may seem quite arbitrary. However, as discussed in Ref. the combination 
of the experimental constraints from the b — > S7 and the B — > tv decays may induce a strong correlation between A 
and mj, restricting their range of variation. This is shown in FigEJ where a scatter plot of these two parameters is 
given and both constraints are applied. Notice that the correlation is further enhanced if configurations plotted in 
grey are considered, corresponding to points falling in region I of FiglH As pointed out in Ref. [2|, constraints from 
rare B-meson decays are affected by uncertainties both in experimental measurements and theoretical estimates, so 
should be considered with care. In particular when the bound from the B — > tv is not implemented the correlation 



1 The quantity Aa M does not depend on the SUSY hadronic sector, so a change in the m<j parameter would imply only a shift of the 
vertical line separating kinematic regions I, III and V from regions IV and VI, without modifying the scatter plot. This implies that for 
rriq> 800 GeV the kinematic region III would be allowed for sequential decays 
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FIG. 5: Scatter plot of the parameters A and raj in the LNM scenario. The color code is the same as in Figure [2] The (red) 
triangle shows the LNM-seq benchmark given in Table ITT1 



benchmark 


Mi[GeV] 


Af 2 [GeV] 


M 3 [GeV] 


ix [GeV] 


tan /? 


m A [GeV] 


mj [GeV] 


m 9 - [GeV] 


mt [GeV] 


A 


LNM-seq 


14 


500 


2000 


126 


34 


97 


300 


700 


444 


-0.08 



TABLE II: The LNM-seq benchmark analyzed in the present paper. 



between A and m t ~ is no longer present [7|. With this caveat, and for the sake of definiteness, we may choose our 
benchmark as a configuration within the grey region of Fig. [5j In particular, in order to kinematically suppress the 
decay Xi ~> ti an d maximize the leptonic sequential signature, we choose as our benchmark a configuration with a 
value of mi close to the upper edge of the allowed range. Our final choice for the last two benchmark parameters is 
plotted in Fig(5]with the triangular symbol, and corresponds to the values given in Table HI1 A=-0.08 and m t -=444 
GeV. 




V. EARLY DISCOVERY OF LIGHT NEUTRALINOS AT THE LHC 



The LHC has already started to put bounds on the supcrsymmetric parameter space. In particular, the very 
constrained SUGRA scenario, in which soft masses and the trilincar coupling are all unified at the GUT scale and 
the fi and ttia parameters are predicted by radiative elecroweak symmetry breaking, appears already to be disfavored 
by the data 4CH44|. Actually, in this scenario lower bounds on both gluino and squark masses are already close to 
the TeV range, in tension with the naturalness picture that is considered one of the motivations of SUSY in the first 
place. 

As already mentioned in the previous section, in order to discuss LHC bounds in the LNM scenario, a dedicated anal- 
ysis was performed in Ref. [?J for some of the specific signatures searched by ATLAS and CMS (namely jets+missing 
transverse energy and one isolated lepton or jets +missing transverse energy and no leptons). As discussed in that 
paper, if squark soft masses of the three families are assumed to be degenerate, the combination of the ensuing LHC 
constraints on squark and gluino masses with the experimental limit on the b — > s + 7 decay imply a lower bound 
on the neutralino mass m x that can reach the value of 11.9 GeV when the gluino mass is at its lower bound, but is 
essentially unchanged for a heavy gluino. However this bound on m x is no longer in place when, as in the present 
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analysis, the universality condition among squark soft parameters is relaxed. This implies that for non-universal 
squark masses the lower bound on the ncutralino mass remains at the value 9 GeV mentioned in the Introduction. 

The LHC is expected to collect C ~5 fb _1 of integrated luminosity at the end of the 2011 run at a center-of-mass 
energy of 7 TeV. In order to estimate the expected signal at the end of the 7 TeV run of the LHC for the LNM- 
seq benchmark introduced in the previous Section we have used ISA JET [46|, applying the same kinematic cuts as 
described in Ref. 40j for the early discovery signature of jets +missing transverse energy and no leptons. The result 
of the simulation is 260 events from SUSY compared to 133 events expected from backgrounds. The background 
estimation is an extrapolation based on CMS measurements. The expected number of SUSY events N is related to 
the SUSY production cross section a = c(pp — > squark) and to the luminosity C by the relation N = ex Cx a, where 
e is the total efficiency due to selection cuts, that for the LNM-seq benchmark we estimate e=0.21. So, in agreement 
to the discussion of the previous Section, the LNM-seq benchmark is expected to provide a slight excess over the 
background, namely at the level of a ~ 3.2 a significance assuming that our estimation on the background has a 5% 
relative uncertainty. 

Such an early hint of SUSY in the 7 TeV run of the LHC would not, however, allow to draw any conclusions on the 
mass and properties of the neutralino, let alone whether the observed excess is compatible to an LNM scenario or not. 
In fact, since the neutralino escapes undetected, its mass can only be reconstructed by observing the sequential decay 
introduced in Eq. ^ where the observed fermions are either muons or electrons in order to have a better discrimination 
of the signal over the hadronic background. In the corresponding final state of 2 jcts+missing transverse energy and 2 
isolated leptons we estimate from the above simulation ~2 signal events at £=5 fb _1 (corresponding to an efficiency 
e=0.00155) obviously insufficient to get any information about the masses. For this kind of analysis E cm =14 TeV 
and a higher collected luminosity will be needed. 



VI. LNM SPECTROSCOPY AT THE LHC WITH £ C a/=14 TEV 

The use of kinematic endpoints to reconstruct the mass spectrum in a sequential decay chain where th e lig htest parti 
cle escapes detection and with the topology as shown in FigQ]has been widely discussed in the literature 36[ H, 39, 47 1 



This technique is based on the simple idea of reconstructing the four unknown masses of the problem {m x , m ; ~, m Xi , rriq) 
by inverting the four observable kinematic endpoints (m^ aa: , rnjlf x , rn^ x o ^ , m ^hi) ) ' wnere 3 indicates the jet produced 
by the emitted quark and mjt(j, ) = Min(mj^,m^,), TUjirhi) = Max(mj.f n ,mjt. ) are defined in order to remove the 
ambiguity between the near and the far leptons £ n , if (shown in FigJTJ) which are not experimentally distinguishable. 
The expected value of such endpoints are given in Table IIIII for the LNM-seq benchmark introduced in Section IIVI 
For completeness we give the analytic expressions of the endpoints as a function of the masses in Appendix [5] 

Two problems however arise in this apparently straightforward procedure. The first issue is related to the fact that 
for particular mass combinations in the sequential decay the four aforementioned endpoints are not independent, since 
the following relation holds [j^: 



K7) 2 = (mTD 2 + K™(«)) 2 - (11) 

In particular, this is true whenever m x < m~/mq. Taking into account the experimental constraint mj > 100 
GeV implies that for a light ncutralino of mass m x = 10 GeV the correlation (fTTj) is verified if rriq < 1 TeV and, in 
particular, holds in the LNM-seq benchmark. To compensate for the occurrence of only three independent variables 
of the physical masses instead of four an additional measurement is needed, for instance the lower kinematic endpoint 
rn Jii{e>iT/^) introduced in Ref. [36| . which corresponds to the lower bound of the rriju histogram with the additional 
constraint (m™ ax ) 2 /2 < (m u ) 2 < {m™ ax ) 2 . 

The second issue related to the analysis of endpoints is that when Eq. (fTTj) is verified the inversion procedure is 
known to have multiple solutions, so that the determination of the masses (m x , m ; ~, m Xi , rriq) is non-unique. This is 
indeed what happens in the LNM-seq benchmark, as shown in Tabic ILTI1 In the top part of the Table we report the 
values for the masses (m x ,m^rri Xi ,rriq). In the middle part of the same Table we give the expected endpoints of the 
kinematic variables (w«, Wj f (; ), m^/^), to^™^^^) in sequential decays, evaluated by using the expressions 
summarized in Appendix [X] By using mass-inversion formulas [39j one finds that starting from the set of the end- 
point values displayed in the Table one recovers the input set of mass values together with a second set of masses. In 
the Table such additional mass spectrum, hereafter referred as the duplicate of LNM-seq, is denoted as LNM-seq'. As 
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Variable 


LNM-seq 


LNM-seq' 


m x 


11 


263 


m l 


305 


383 


m Xi 


515 


688 


rriq 


703 


896 


max 
III go 


415 (417.5± 3.5) 


m max 


632 (631.3± 3.8) 


max 
m jl(lo) 


338 (342.2 ±4.3) 


m max 


477 (483± 14) 


min 
' n jU(e>ir/2) 


400 (399.3 ± 1.7) 


n',p> 


282,385 


232,477 



TABLE III: Expected endpoints (in GeV) of the kinematic variables (mee,mjee,mji(i ),mj e ( hi ),m'J}l™ e>7T /2)) f° r the LNM-seq 
benchmark of Table [Tit calculated using the expressions summarized in Appendix [X] In parenthesis we give measurements of 
the same quantities analyzing the output of a simulation of proton-proton collisions at Ecm = 14 TeV assuming an integrated 
luminosity of 100 fb _1 (see text). The benchmark LNM-seq' indicates the duplicate model of LNM-seq, i.e. a different mass 
pattern providing the same values of the observed endpoints in one- dimensional distributions [isj. In the last line the coordinates 
of the point (m,je(i ),mj£(hi))=(n' ,p') can break the degeneracy between the two duplicated models. For this particular mass 
pattern it is not possible to measure the exact values n' ,p' , since they lie on a straight boundary of the two-dimensional 
distribution (see Fig. 9(a) I. However this is sufficient to break the degeneracy (see text). For this reason we do not provide a 
measured value of n! ,p' from the simulation. The analytic expressions of (n',p') as a function of the physical masses are given 
in Appendix IA1 



will be discussed in the following, in order to overcome this duplication problem it is necessary to go beyond endpoints 
in one-dimensional histograms and to analyze the correlations among different invariant masses in two-dimensional 
plots 0. 

We wish now to discuss if the procedure outlined above can be applicable to determine the mass spectrum of the 
LNM-seq benchmark using the LHC data at Ecm=^ TeV. In order to do this, we simulate proton-proton collisions 
at i?c7M=14 TeV using ISAJET and select events with two jets, two isolated leptons and missing transverse energy. 

Notice that fast detector simulation tools which have been developed for the study of specific Supcrsymmetric 
scenarios such as in Supergravity-inspired (SUGRA) benchmarks are not available for the model under consideration 
here. So, in order to take into account the detector response, and specifically the uncertainty in the reconstruction of 
jet energies, in our simulation we smear the energy E of quarks and gluons in the final state. In particular, we apply 
a resolution which depends on energy as 0.9 \[E. The Missing Transverse Energy (MET) is then determined by the 
vector sum of the energies of the neutrinos and the LSP plus any smearing applied to the hadronic jets. Due to the 
large mass of squarks and gluinos the dominant background is expected to be due to tt production. In particular, 
as shown below, large cut values are needed for an effective separation between the signal and the background. We 
apply to the output of our simulation the following cuts: 

• The two leptons (e + e~ or are required to satisfy |?7| < 2.4 (where r\ = — ln[tan(6*/2)] is the pseudo- 
rapidity and 9 is the angle with the beam axis) and pr > 20 GeV (where px is the transverse momentum) . 

• The kinematic separation between outgoing states is required to be AR > 0.5, where AR = {Arf) 2 + (A<fr) 2 
and <j) is the azimuthal angle. 

• A missing transverse energy Et > 300 GeV is required, in order to indicate the presence of high energy 
neutralinos. 

• The scalar sum of the transverse momenta pr of leptons and jets is required to be larger than 600 GeV. 

• In the study of 2-dimensional distributions we have removed events where the invariant mass of the two outgoing 
leptons falls in the range 87 GeV < Mu < 97 GeV, in order to subtract Z boson decays. These events include 
also those produced in the branched decays shown in FigJTJ right. 
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FIG. 6: Event distribution as a function of the final leptons invariant mass mu for the benchmark LNM-seq. Plots on the 
left show the dilepton invariant mass at the preselection level while plots on the right show the same histograms after the cuts 
described in Section fVTl are applied, with the exception of the subtraction of the Z peak (which is clearly visible). Upper plots 
show the case when the final leptons are given by e + e~ and fi + p~ , while lower plots show the r + r~ final state. In the plot on 
the upper left the hatched histogram shows the distribution of the ee + up events, while the (green) solid histogram shows the 
same for the ep + pe events, which provides an estimation of the expected contribution from the SUSY background (see text). 
The shaded histogram in the plot on the upper right shows the background-subtracted distribution given by the difference 
between the ee + fip and the ep + /ie histograms. In the case of the lower plots the hatched histogram on the left shows the 
distribution of t + t~ events before cuts, while that on the right shows the same quantity after cuts. In all plots the white 
histogram shows the ti backgrounds. 



Events with more than two jets or two leptons are rejected to minimize the effects of combinatorics. Out of the 
two jets in the event, only one must be associated to the dilepton in order to construct the rrijm \, m^n^hi) invariant 
masses. The jct-dilcpton pairing is found by choosing the combination that yields the smallest value for the muj 
invariant mass. 

The expected distribution of the events with 2 jets+2 leptons+missing transverse energy produced in the simulation 
is shown in Fig(5] as a function of the invariant mass m« and in Fig(7] as a function of muj ■ In both Figures the plots 
on the left show events at the pre-selection level, while those on the right show the same distributions after the cuts 
described above, with the exception of the subtraction of the Z peak. Moreover, upper plots show the case when the 
final leptons are electrons or muons, while lower plots show the t + t~ final state. In all plots the white histogram 
shows the ti backgrounds, which, as can be seen, is strongly suppressed by the cuts. 

A different and potentially sizeable source of background making the determination of endpoints difficult is also 
represented by SUSY events where the two charged leptons used to calculate the invariant mass are not originated 
in a sequential decay, but are produced instead by the decays of charginos originating from different decay chains. 
In the upper-right plots of Figs. |B] and [7] these undesired events are subtracted exploiting the fact that in this case 
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FIG. 7: Event distribution as a function of the invariant mass meej of the final leptons and the jet for the benchmark LNM-seq. 
The color code and the cuts are the same as in Fig(5] 



the flavors of the two leptons is uncorrelated, while when the two leptons are produced in the same sequential decay 
they have the same flavor. For this reason an effective subtraction of this background is obtained by taking the 
difference between the number of ee+fip same-flavor events minus the number of events where the flavor of the final 
leptons is different, ep+fie. This is indeed an effective technique to subtract chargino decays and to allow a better 
identification of the endpoints, since the difference of the two distributions is expected to drop beyond the boundaries 
of the sequential process. In particular, in the plots on the upper left of Figs[S]and[7]the hatched histogram shows 
the distribution of the ee + fip events while the solid histogram shows the same for the ep + fie events, which provides 
an estimation of the expected contribution from the SUSY background. Moreover, the shaded histogram in the plots 
on the upper right show the background-subtracted distribution given by the difference between the ee + fip and 
the ep + fie histograms. In the case of the lower plots the hatched histogram on the left shows the distribution of 
t + t~ events before cuts, while that on the right the same quantity after cuts. Comparison of left-hand figures and 
right-hand ones prove the overall effectiveness of the applied cuts to subtract the standard model background, in 
particular from top decays. 

Notice that Figs. [5] and [7] are normalized to the luminosity and represent the theoretical expectations of the 
corresponding distributions. They are obtained using an integrated luminosity of 546 fb^ 1 (or 3940 events after 
selection), allowing in particular to determine easily the position of the endpoints (that agree with the values given 
in Table IIIip . However, for a lower value of the integrated luminosity the position of the endpoints is blurred 
by statistical fluctuations, worsening their determination. In the following we will assume for an optimistic and yet 
realistic prediction of the latter quantity £=100 fb _1 at the end of the 14 TeV LHC run. The corresponding prediction 
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for the mu distribution for a simulated experiment is given in Fig. 8(a) where the ee+/i/x-e/x — fie subtraction between 



same-flavor and different flavor events has been applied to reduce the SUSY background from chargino decays. After 
selection cuts but without subtracting the Z peak this plot contains 726 events, which become 502 when the Z peak 
is subtracted. 

The endpoint of this histogram provides the first edge m 1 ^ ax needed for the kinematic reconstruction of the masses. 
In order to find it we employ a method inspired by an edge-detection algorithm frequently used in the field of image 
processing and computer vision 48| . Actually, the most sensitive method for measuring the position of an edge would 
be to obtain its expected distribution from a simulation and to perform a likelihood fit to the data. However, this 
method depends on SUSY parameters, hence lacking generality. On the other hand the edge-detection algorithm, 
while it may not yield the best sensitivity, has sufficient generality to be applied to a wider range of problems. 

The edge-detection algorithm is a method that allows to find the endpoint of a sharply falling distribution by filtering 
the data histogram through an appropriate function. For concreteness, we take as a filtering function f(x,fx,a) = 
2sinh((x — fi)/a) / cosh 3 {{x — n)/a) and try to minimize the quantity: 



N da 



(12) 



with respect to /z and with a fixed. In Eq. (jT2"|) Xi represents the data count in the i-th bin of the histogram. The 
width parameter a has the effect of smoothing the distribution, hence making the algorithm immune to noise. The 
choice of a is determined by looking at the width of the distribution. If the value is too large, the edge determination 
is imprecise, whereas if the value is too small, then it will be sensitive to outliers. 

In Fig. |8(b)| we apply the method outlined above to filter the m« histogram of Fig 8(a) We assume cr=30 GeV and 
the endpoint position is represented by the rightmost minimum, and is reported on the data histogram in 8(a) with a 
vertical solid line. In order to estimate the statistical fluctuation of m™^ x we then repeat the same procedure for 100 
pseudo-experiments identical to the one analyzed in Figs. 8(a) and |8(b)"| The corresponding frequency histogram for 
the outcome of m"\ ax is given in Fig 8(c) In this way we find TOj,™ a:c =417.5± 3.5 (this value is reported in parenthesis 
in Table |nl|. 

In order to find the other endpoints needed to reconstruct the masses, in principle the above procedure can be applied 
also to the histograms obtained by plotting the same simulated events as a function of the other invariant masses muj , 
mjt(hi) and mji(io)- However, in the latter distributions the position of the endpoints cannot be determined accurately 
because the number of events is not large enough to saturate the endpoint of the histogram, which systematically 
drops at a value considerably lower than the true one for a lack of points in the tail. In this case an unambiguous 
determination of the endpoint is strictly speaking impossible, and only some educated guess can be made. In order 
to do this it can be useful to resort to two-dimensional plots. This is done in Figs 9(a) and 10(a)| where the events 
of the pseudo-experiment plotted in Fig 8(a) and that lie to the left of the determined value of m^ ax , are plotted in 



the planes m^a^-m^r^and mu-rriju. In this way both plots contain 497 events. 

In the case of the TUjia^-rrij^fh^ two-dimensional plot of Fig. 9(a) the shape of the region covered by the data 
points nicely fits an isosceles triangle. This very symmetric shape is expected in two situations [39| : i) if the slepton is 
produced off-shell (i.e. if mj > m Xi ), since in that case there is no longer distinction between the near and far lepton 
and the two leptons have exactly the same kinematic properties; ii) when the following relation among masses holds: 



2 
I 



< 



< 



> / 2 



(13) 



and at the same time there is a large hierarchy between the slepton mass and the neutralino, m x /mr <Cl. In the latter 
case, that corresponds to our LNM-seq benchmark, the expected boundary is actually delimited by four vertexes, but 
the fourth point (n',p'), that is supposed to be used to break the degeneracy among duplicate models and whose 
coordinates are given in Table UTTT lies on the straight line n! + p ' =(mju max ) 2 when m x <C mj (it is represented by one 
of the two filled circles in Fig 9(a) ) and cannot be observed 2 . Notice however that, in spite of this, the degeneracy 



2 The density of points of the two— dimensional distribution is expected to have a step-like drop for Wj£(; )> p' [39ll , allowing in principle 
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FIG. 8: (a) Histogram of events with two jets, two isolated leptons and missing energy after cuts (but including the Z peak) 
as a function of mu for the LNM-seq benchmark in a simulated experiment at the LHC with y/s=14 TeV and assuming 100 
fb _1 of integrated luminosity. The plot contains 726 events (502 applying the cut on the Z peak), (b) Data filtering of the 
data in (a) through the function given in Eq. (|12[) . plotted as a function of the edge position guess and for a fixed value of 
the width parameter a = 30 GeV. The endpoint m^} ax of the distribution is represented by the rightmost minimum, and is 
reported on the data histogram in (a) with a vertical solid line, (c) Frequency histogram for the outcome of rrif} ax for 100 
random pseudo-experiments identical to the particular one plotted in (a) and (b). 



a determination of p' . We have verified that in practice this measurement is not possible because of the large fluctuations in the 
determination of the density due to the low number of points. 
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FIG. 9: (a) Two-dimensional distribution in the plane 
left of the value of mS" 



determined in Fig. 8(b) 



+li)j(lo) m (e+ M )j(hi) 



of the events plotted in Fig 8(a) that lie to the 



and when the cut on the Z peak is applied. The plot contains 497 events. 
The (red) solid line represents the expected boundary for LNM-seq, while the (blue) dashed lines are the fits of the boundaries, 
when the shape is assumed to be a triangular one (see text). The shaded area is the expected boundary for the duplicated 
model LNM-seq'. The two filled circles show the positions of the point (n',p') for LNM-seq and LNM-seq'. Although for 
LNM-seq the exact position of (n',p') cannot be measured, the fact that it lies on the boundary of the triangular shape is 
sufficient to break the degeneracy between the two duplicate mass patterns LNM-seq and LNM-seq' (see Table [TTTI and text), 
(b) Frequency histogram for the output of the quantities rnj^^ and rn™^^ for 100 random pseudo-experiments identical to 
the particular one plotted in (a). In each pseudo-experiment mjl^ and rn™^^ are obtained as the crossings of the fitted lines 
with the axes. 
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with the LNM-seq' model is easily broken, since the expected shape for the LNM-seq benchmark, represented in Fig. 



9(a) by the (red) solid triangle, is very different from the corresponding one for LNM-seq' shown as the shaded area 



in the same figure. As a consequence of this the following relation among the masses: 



< — i- < _2l 



in v 



in r 



(14) 



which corresponds to the particular trapezoidal shape of the LNM-seq' benchmark, can be safely discarded. 

The distribution of points in Fig 9(a) is clearly not dense enough to saturate the vertexes of the triangle 3 , so strictly 



speaking other kinematic regions different than the LNM-seq benchmark cannot be ruled out (for a summary of shapes 
corresponding to different kinematic situations see for instance Fig. 8 of [1^]). However the shape is very compatible 
to a isosceles triangle as in the LNM-seq benchmark. Moreover in the LNM-seq case the possibility that the triangular 
shape is due to an off-shell sequential decay can be easily excluded on dynamical grounds. In fact the branching ratio 
of the off-shell sequential decay drops by at least two orders of magnitude compared to the on-shell situation. In this 
case, in order to detect a few hundreds events in the sequential channel as in Figs 9(a)| 10(a) the production cross 

This in turns would 
In the LNM-seq 



section would need to be much larger than in the LNM-seq, say in the range of a few tens pb. 
lead to a dramatic enhancement of branched decays that in this case would be the dominant ones, 
scenario this would lead to a huge number of events showing up in the Z or Higgs peaks when plotted as a function 
of m«. The non observation of such an excess would easily allow to rule out that the events plotted in Fig 9(a) are 
due to off-shell decays. 

In such a predicament we then propose to make the educated guess that the shape is a triangle and that the 
sequential decay is on shell. In this case, assuming that the boundaries of the region are straight lines, the edge- 
detection method that we used for the mgg one-dimensional histogram can be modified to find the position of the 
edges in the two-dimensional x-y plane by minimizing: 



± 
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(15) 



The sign '±' should be chosen depending on the observed slope of the boundary. 

The result of the above procedure is shown in Fig. 9(a)| where the fitted straight boundaries are represented by the 
dashed (blue) lines. In this figure the two corresponding endpoints m^n^y ml Ji(hi) are then obtained as the crossings 
of the boundaries with the two axes. In Fig {9(b)1 thc same procedure is repeated for 100 pseudo-experiments identical 

md to^;^) are given. The 

is reported in parenthesis in Table ITlTl 

is finally discussed in Fig. [TDJ As in the previous 



to the one shown in Fig 9(a) and the frequency histogram for the output values of m™ 1 ^ 
corresponding determination for m™?£i and m mc 



max 

je(io) " l j£(hi 

The determination of the last two endpoints mj^ x and ''"-^^(e <7I ./2) 
figure the (red) solid line represents the expected boundary for the LNM-seq benchmark, while the shaded area is the 
corresponding one for the duplicate model LNM-seq'. Also in this case the upper and lower boundaries of the region 
covered by the simulated data are compatible with strai ght lines. This is broadly consistent with the guess that the 
kinematic region of Fig 10(a) is due to an on-shell decay (39j. When the boundaries are intersected with the value of 



in/ 



determined in Fig 8(b) 



hnth m max nnH m mm 

Dotn m ]U ana m jU ^ e< ^ /2) 



can be obtained. Finally, in Fig |10(b)| repeating the same 



procedure for 100 pseudo-experiments identical to that of Fig 10(a) the frequency histograms for the output values 



m 1 : 



m 



(9<7r/2)> 
1 min 

"JU — ■' l j«(6»<7r/2) 

Once the 5 endpoints m 



and of the slopes of the upper and lower boundaries is obtained. The ensuing determinations 
are given in parenthesis in Table Mil 

are determined (notice that, as mentioned before, 



<7 r/2: 



since the relation of Eq. (fTTj) holds, only four of them arc independent) they can be used to determine the masses. The 
mass inversion is obtained in a straightforward way by simulating a large number of random values of the four masses 
m Xi , TJiq and plotting the histogram of the mass combinations whose theoretical values of the endpoints 



/■/? 



v 



fall within the measured ranges. The result of such an inversion is shown in Fig. 1111 As expected, this procedure 



If it were so the endpoints would be observable in the correspondent one-dimensional projections of the distribution, without the need 
to resort to two-dimensions in the first place! 
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The (red) solid 



FIG. 10: (a) Two-dimensional distribution of the same events of Fig. 
line represents the expected boundary. The (blue) dashed lines are the fits of the upper and lower boundaries when they are 
assumed to be straight lines (see text). The value of m^} ax is the one determined in Fig. 8(b) The shaded area is the boundary 



for the duplicated model LNM-seq'. (b) Frequency histogram for the outcome of the slopes of the upper and lower boundaries 



of the region in (a) and of the endpoints m!^ x and m£ 



9>7r/2) 



, for 100 random pseudo-experiments identical to the particular 



one plotted in (a) . The latter quantities are obtained as the crossing points of the relevant fitted boundary lines. 
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Mass determination for the particles of the sequential decay obtained from the values of the endpoints measured 



FIG. 11 

using the simulated sample of Figs 



9(a) and 10(a) Upper left 



m. 



upper right: 



m,~; lower left: 



It!, 



lower right -.nig. In 

all plots the histogram bins are filled with a shaded box in correspondence to mass patterns verifying Eq. (|13[) . Only the 
mass determinations corresponding to the shaded peaks are compatible to the data, namely to the shape of the data region of 
Fig|9(a)|in the m je{hi y 



ojplane. 



leads to two different solutions, corresponding to the LNM-seq benchmark and to the duplicate one LNM-seq'. Notice 
however that from the discussion of Fig 9(a) this degeneracy can be easily broken. In fact the distribution of the 
simulated data points in the TO^^-m^^i) plane is strongly inconsistent with the LNM-seq' solution, allowing to 
conclude that only mass patterns verifying Eq. (|13[) are compatible with the simulated data. In Fig llll such mass 
patterns are marked by filling the bins with a shaded box. From this figure one can see that the correct solution can 
be clearly discriminated from the duplicate one for all the masses involved the decay. In this way from Fig. QT]we get 
the following determination of the masses: 



(103 ± 43) GeV 
(349 ± 27) GeV 
(561 ± 28) GeV 
(751 ±28) GeV. 



(16) 
(17) 
(18) 
(19) 



The reconstructed value for the ncutralino mass deviates remarkably from its input value. This is explained by 
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the fact that the center-of-mass energy available in the sequential decays is set by the squark and is much larger 
than the ncutralino mass. So neutralinos are produced in the relativistic regime, in which their kinematics is almost 
insensitive to the actual value of m x . Moreover, the dependence of the kinematic endpoints on the neutralino mass 
is quadratic, so determining m x implies taking a square root, which worsens the accuracy compared to that with 
which the endpoints are measured. However, the reconstructed value for m x , though deviated from the actual value, 
would entail the important indication that some neutral stable particle is being produced in the process, whose mass 
can be compatible to a light range. Notice that the similar yield to leptons of the three families observed in Figs. [5] 
and [7] allows to conclude, as discussed in Section 1111 Al that both the neutralino and the Xi particle are of gaugino 
type. As explained in Section [IT] this implies m x ~ Mi and m Xi ~ M2, so the reconstructed mass spectrum would be 
incompatible with neutralinos in SUGRA scenarios, since it would point toward a ratio Mi/M 2 < 0.5, in disagreement 
with what expected in models where gaugino masses are unified at the GUT scale. 



VII. CONCLUSIONS 



Signals of relic particles in direct DM searches raise the interest for masses of these relic particles in the range 7-8 
GeV < m x < 50 GeV. This is actually the mass range compatible with the annual-modulation effect measured by the 



DAMA collaboration ll| , when this effect is interpreted in terms of DM particles with an elastic coherent interaction 
with nuclei. If the positive results of other experiments of DM direct detection are taken into account (CDMS 0|, 
CoGeNT 0, CRESST 0) the mass range restricts to 7-8 GeV < m x < 15-20 GeV Q. These experimental results 
arc fitted quite well by the Light Neutralino Model [H, Q that is an effective Minimal Supersymmetric extension of 
the Standard Model at the electrowcak scale without requirement of a gaugino-mass unification at a grand unification 
scale. At variance with Supergravity-inspired (SUGRA) models, in the LNM the neutralino mass can be as small as 
about 9 GeV, as discussed in the Introduction. 

In view of the interest of these light neutralinos in the phenomenology related to DM direct detection, in the present 
paper we have addressed the problem of a search at the LHC for a neutralino of very light mass. A preliminary analysis 
in this direction was performed in Ref. [I3|. There, specific scenarios and benchmarks within the LNM and dictated 
by relevant cosmological properties were considered, and the relevant expected event rates determined; however, no 
specific analyzes of the signal/background ratios and of kinematical distributions were performed. 

In the present paper we have extended the investigation of Ref. [lit by making use of numerical simulations to 
estimate in a realistic way the detectability of light neutralinos at the LHC over the SM background and to show what 
information about the masses of SUSY particles can be extracted from the data. Within the sequential and branched 
decay chains that constitute the typical processes by which neutralinos can be searched for at the LHC, we have 
singled out the most dominant contributions in the context of the LNM. For this aspect, the analytical expressions 
for the light neutralino spectroscopy reported in Sect. II have been exploited. 

We have then selected a benchmark (dubbed LNM-seq) with the specific feature of belonging to a value (around 
10 GeV) of the ncutralino mass that is representative of the neutralino mass values in the low-side of the m x ranges 
mentioned above. In terms of this benchmark we have investigated which are the expectations for having some signals 
at the LHC at different stages of the LHC operation. 

We have found that with the integrated luminosity C ~5 fb _1 that would be collected at the end of the 2011 run at 
a center-of-mass energy of 7 TeV, the LNM-seq benchmark is expected to provide a slight excess over the background, 
namely at the level of a ~ 3.2 a significance, assuming that our estimation on the background has a 5% relative 
uncertainty. As discussed in Sect. [V] this would not be sufficient to draw any conclusions on the mass and properties 
of the neutralino. 

We have then analyzed the prospects in terms of the integrated luminosity of C ~100 fb _1 that might be reached 
by the LHC at the end of its 14 TeV run. To this purpose a detailed analysis has been performed by employing 
one-dimensional and two-dimensional distributions with the scope of establishing how the inputs of the LNM-seq 
benchmark can be reconstructed by the determinations of the relevant end-points in the various mass distributions. 
The problem of the disentanglement of the true solution from the duplicate solution in the inversion procedure from 
the end-point values to the model parameters has been addressed. 

The main result concerns the reconstruction of the neutralino mass that finally turns out to be determined as 
(m x )rec — 103 ±43 GeV. This value deviates remarkably from the input value for m x , but this result is not surprising 
in view of the difficulty in reconstructing the mass of light stable particles in relativistic events. However, the 
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reconstructed value for m x , though deviated from the actual value, would entail the important indication that some 
neutral stable particle is being produced in the process, whose mass can be compatible to a light range. Moreover, 
the reconstructed masses would suggest Mi/M 2 ^0.5, in disagreement with what expected in SUGRA models where 
gaugino masses are unified at the GUT scale. 
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Appendix A: Analytic expressions of kinematic endpoints 

We give here for completeness the analytic expressions used to calculate the kinematic endpoints in Tabl dllll These 
formulae are taken from Ref.l39j|. 



(mZ ax ? 
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with: 
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(mJC) 2 = m 2 9 (l-R Xi5 )(l-R lxi ), 
(mf e a f x Y = rof (1 - R Xi9 ) (1 - R xI ) , 

(Kit^y = m \ a-^ja -^r) ( 2 - V 



(A6) 
(A7) 
(A8) 



and i?; m = mf/m^ with i = x, i, Xii 9- 

Moreover, the quantities (n',p') are given by: 



min(n,p) 
max(n,p), 



(A9) 



where: 



m?(l-i? Xl9 ~) 



V/ 



(A10) 
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